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We study with exact diagonalization the zero temperature properties of the quarter-filled extended 
Hubbard model on a square lattice. We find that increasing the ratio of the intersite Coulomb 
repulsion, V, to the band width drives the system from a metal to a charge ordered insulator. The 
evolution of the optical conductivity spectrum with increasing V is compared to the observed optical 
conductivity of several layered molecular crystals with the 9 and /3" crystal structures. 
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Charge ordering in strongly correlated electron systems 
is currently under intense investigation. Charge order- 
ing is relevant to a broad range of materials including 
the cuprates[[l|, manganatesQ, magnetite ||, vanadium 
oxides [Q, and the Bechgaard salts ||. The 9 and f3" types 
of layered molecular crystals based on molecules such as 
BEDT-TTF [= bisethylenedithio-tetrathiafulvalene] §, 
display charge ordering, metallic, and superconducting 
phases close to each other JjJ. Charge ordering driven 
by a strong inter-site Coulomb repulsion || || is possible 
in crystals with the 9 and (3" arrangements of BEDT- 
TTF molecules because their bands are quarter-filled 
with holes, in contrast to the well studied K-type, for 
which strong dimerization of the molecules lead to a half- 
filled band . The #-type crystals undergo a transition 
from a metal to a charge ordered insulator as the temper- 
ature, pressure, uniaxial stress, or anion is varied [Q, ||. 
Furthermore, the metallic phase exhibits features charac- 
teristic of a strongly correlated system. In particular, the 
optical conductivity spectra display a broad mid-infrared 
band and a near absence of a Drude-like peak|ll|. This 
is in contrast to conventional metals, for which the total 
spectral weight is dominated by a Drude peak. 

In this Letter, we use the results of an exact diagonal- 
ization study of the relevant extended Hubbard model to 
argue that the inter-site Coulomb repulsion is responsible 
for the observed metal-insulator transition in the 9 and 
(3" crystals. We show how the Drude weight decreases as 
the inter-site Coulomb repulsion, V, is increased, until, 
at a finite value of V, a transition to an insulating phase 
occurs. Simultaneously, long range charge ordering grad- 
ually sets in. We further find that a redistribution of the 
optical conductivity spectra occurs close to the metal- 
insulator transition. This finding is in qualitative agree- 
ment with experimental data on 9 and (3" organic salts 

The quarter-filled extended Hubbard model on a 
square lattice is the simplest strongly correlated model 
that can potentially describe the competition between 



metallic, superconducting, and insulating phases in the 
9 and (3" materials . The Hamiltonian is 

H= -t ]T (4<7 C 3<r + c] a C la ) + n A n ii 
< ij > , cr i 

<ij> 



where c ia creates an electron of spin a at site i. For 
V = 0, previous calculations suggest that the system 
is metallic with no charge order ]l7j]. In the limit of 
f/,V>t the double occupation of sites is suppressed and 
the ground state is insulating with checkerboard charge 
ordering and long range antiferromagnetic correlations 
along the diagonals . 

In the present work we consider large (but finite) U /t 
and vary V/t fll8[| . Previously, an SU(N) generalisation 
of the model with U — > oo, was studied in the large 
N limit using slave bosons jl5| It was found that 

as V/t increased there was a transition from a metallic 
phase to a superconducting phase (with d xy symmetry) 
to a charge ordered phase. Given this potentially rich 
phase diagram it is important to determine whether or 
not these previous results are an artefact of the approx- 
imations used or whether they reflect the actual physics 
(N=2 and finite U/t). Exact diagonalization of small 
systems provides such a test. Previously, similar models 
with more parameters, aimed at a more realistic descrip- 
tion of the details of the materials, have been studied 
by Hartree-FockQ, and quantum Monte Carlo[[j"9|, and 
exact diagonalization |?0| . We performed a Lanczos cal- 
culation on L = 8, 16, 20 site clusters[[l7j avoiding un- 
compensated spin moments in the cluster (jVj = iVj_ ) 
at quarter-filling: < n >= 1/2. A powerful method to 
determine from small systems whether the bulk system 
is metallic or insulating is to evaluate the Drude weight 



2 




1 f 2 3 4 5 



V co /t 

c 

FIG. 1: Metal-insulator transition induced by the nearest- 
neighbour Coulomb interaction V. The Drude weight, D, is 
shown as a function of V/t, for L — 16 and various values of 
U/t (continous and dotted lines) and for L = 20, U = lOt 
(dashed line). The arrow in the horizontal axis marks the 
onset of checkerboard charge ordering at V = V/?° , for U = 
lOt (see Fig. ^|). The inset shows the finite-size scaling of D 
as a function of 1/ L x for different values of V/t with U = Wt. 
The metal-insulator transition occurs at V C MI « 2.2t. 

D HI, H ||. It is given by 

J _ {0jT|0) 1 v |(n| J:c |0)| 2 
27re 2 4L L^E n -E { ' 

where i?o and _E„, denote the ground state and excited 
state energies of the system, respectively. T is the kinetic 
energy operator and j x is the current operator in the x 
direction at zero wavevector (q = 0). The occurrence of 
an insulating phase is marked by the exponential van- 
ishing of D with the linear size of the system L x = \f~L 

In Fig^jl] we plot the Drude weight as a function of V/t, 
for L = 16 and different values of U/t. For U = lOt we 
also plot the Drude weight for the cluster with L = 20. 
As the Coulomb nearest-neighbors repulsion is gradually 
increased, the Drude weight decreases until it eventually 
vanishes. For L = 16 we estimate the critical value for 
the metal-insulator transition to be ~ 2.2i; for this 
value, the Drude weight becomes of the order of 10~ 3 
and stays of this order up to V = 5t. For L = 20 we do 
not find any significant change in the dependence of the 
Drude weight with V, meaning that the finite size effects 
in the Drude weight are weak. This is confirmed in the 
inset of Fig. [I] which shows the finite-size scaling of the 
Drude weight. Indeed, we find that for V > V C MI « 2.2t, 
the Drude weight displays an exponential dependence 
with l/L x , as expected for an insulator(22j |3j. It ex- 
trapolates to values of the order of 10~ 5 in the thermo- 
dynamic limit. In contrast, for V < V c , the Drude 
weight is weakly dependent on 1/L X extrapolating to a 



finite value in the thermodynamic limit, consistent with 
a metallic state. 

The occurrence of charge ordering can be investigated 
by computing the charge correlation function 

C (V) = lJL e iq(R '- Rj) < n inj > (3) 

ij 

where q are the allowed momenta on the cluster. We find 
that at V — 0, C(q) is featureless. As V is increased, 
C(q) peaks at q = Q = (n, it), signalling checkerboard 
charge ordering. In our calculations we do not find any 
evidence for other charge ordering patterns, such as those 
found in References |], ^0|. The system has long range 
charge ordering if C(Q) remains finite in the thermody- 
namic limit. In Fig. g we plot C(Q) as a function of l/L. 
A linear dependence with l/L is expected for the finite 
size scaling of an order parameter with a discrete sym- 
metry in two dimensions [p5| . For V = 0, we observe that 
the linear extrapolation ofC(Q) (using L = 16, 20) to the 
infinite volume limit tends to zero, i.e., the charge is ho- 
mogeneously distributed in the lattice. We find that the 
extrapolation of a linear fit of L = 8, 16, 20 to the ther- 
modynamic limit becomes finite at about V^° w 1.6t 
p6[ . If we decrease U the critical value V^° increases. 
This is because doubly-occupied sites proliferate so that 
inducing the checkerboard charge ordered state becomes 
energetically less favorable. Further insight into the tran- 
sition can be obtained by computing the charge order- 
ing parameter [p7|, rj = y / C(Q)/ < n 2 >. A non-zero 
•q means long range charge ordering. In particular, for 
X] = 1 the checkerboard charge ordered state is fully 
formed, while for r\ — the charge is homogeneously 
distributed in the lattice. At V — 2.2t we find that the 
checkerboard is close to being completed, 77 w 0.75. We 
note that the results we find here are different from the 
ones obtained in two-leg ladders |27| where a transition 
from an homogeneous insulating phase to a charge or- 
dered insulating phase takes place. This difference is due 
to the one-dimensional nature of the two-leg ladder. 

At this stage we can provide a physical interpreta- 
tion of our findings. Only when V is sufficiently large 
so that the checkerboard charge ordered state is nearly 
completed, a transition to an insulating state is possible. 
This is because moving an electron within the checker- 
board ordered state would cost an energy of « ZV , which 
is comparable to the bandwidth, W — 8t. As a conse- 
quence electrons can become localized. As an aside we 
note that Fig. |l| shows that increasing U from lOt to 20t, 
leads to only a slight decrease in the critical value, V^ n . 
This contradicts the results obtained by Ohta et al. |28[ ] 
who claim that V c « J7/4. We note that the critical value 
yMi £ or i ar g e jj follows the tendency encountered in one- 
dimensional rings ^| and two-leg ladders [^7) for which 
V C M/ decreases to 2t as U -^00. 

Since we find that V C M/ > V^° , it might be pos- 
sible that a charge ordered metallic state is realized 
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FIG. 2: Checkerboard charge ordering induced by the near- 
est neighbour Coulomb interaction. Finite-size scaling of the 
charge correlation function C(q) (defined in Eqn. (^)) at 
q = Q = (tt, 7t) for U = Wt and several values of V/t at quar- 
ter filling. Long range checkerboard charge ordering occurs 
for V > Vf° » 1M. 



in the range: Vf° < V < V C M/ . However, cau- 
tion is in order because the former value is quite sen- 
sitive to finite-size effects. A possible candidate for this 
phase is the quarter- filled organic crystal: /3"-(BEDT- 
TTF) 2 SF 5 CH 2 CF 2 S03, which displays metallic behavior 
of the resistivity and at the same time charge dispropor- 
tion in alternate molecules [[30|. 

In order to further our understanding of the metal- 
insulator transition and make contact with experiments 
on 9 and 0" materials, we have also computed the real 
part, cr(u>), of the optical conductivity (3^] at frequency 



a(cu) = DS(lo) 



ire 

T 



E n — Eq 



E n + Eq). 



It obeys the following sum ruleBS 



f°° ire 2 

/ a(uj)duj = -— < 0|T|0 > 

Jo 



(4) 



(5) 



In Fig. U we plot the evolution of the optical con- 
ductivity for different values of V/t and L = 20. We 
also found that the plots were qualitatively similar for 
L = 16 (not shown). For the regular Hubbard model, i.e. 
V = (Fig. U (a)), there are two features at non-zero 
frequency. The highest energy feature is related to tran- 
sitions associated with the Hubbard bands and the other 
is the "mid-infrared" band, previously found by Dagotto 
et aZ.@. At V = IM (Fig. 1(b)), the Drude peak is 
reduced (see Fig. [l]), signalling the incipient localization 
of the charge carriers. Spectral weight is transferred from 
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FIG. 3: Transfer of spectral weight in the optical conductivity 
with increasing V/t. The frequency dependence of the real 
part of the conductivity (t(lu) for L = 20 and U = lOt is 
shown in both the metallic [V = (a), V = 1.8t (b), V = 2t 
(c)], and insulating [V = 3t (d)] phases. Note that close to 
the metal-insulator transition the spectral weight from the 
Drude peak and the Hubbard feature is transferred to the 
"mid-infrared" peak. 



both the Hubbard feature and the Drude peak to the mid- 
infrared band found at V — 0. This effect is more clearly 
observed for V = 2t in Fig. || (c): the Drude peak is fur- 
ther suppressed and most of the spectral weight appears 
at low frequencies, in the range V < lo < 3V. Increasing 
V further to just below Vmi gives a similar shape of the 
optical conductivity (not shown) with a very small Drude 
weight. Finally, Fig. || (d) shows the optical conductiv- 
ity on the insulating side of the transition: V = 3i. An 
optical gap opens and most of the spectral weight con- 
centrates at larger frequencies building up a single broad 
resonance. If V is further increased, V >> U/4, then 
the band associated with charge excitations due to V ap- 
pears above the Hubbard band situated near u) ~ U, as 
discussed by Ohta et al. |^8| . In summary Fig. 3 shows 
that close to the metal-insulator transition the optical 
conductivity is dominated by the "mid-infrared" feature, 
for both metallic and insulating phases. We now discuss 
how this is what is observed in the 9 and (3" materials. 
The optical conductivity we find close to the 



metal-insulator transition, V ^ V, 



MI 



is qualita- 



tively similar to the optical conductivity of /3"-(BEDT- 
TTF) 2 SF 5 CH 2 CF2S03, measured by Dong et al. 
(at a temperature, T = 14 K). Indeed, the small 
spectral weight at zero frequency found here is consis- 
tent with their failure to observe a Drude peak, even 
though the system is metallic. Also their observation 
of large spectral weight at low frequencies is consis- 
tent with the mid-infrared band we find in our calcu- 
lations. At V < V^ n , the system is still on the metal- 
lic side of the metal-insulator transition which is con- 
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sistent with the observed metallic behavior of the re- 
sistivity. From the above discussion we conclude that 
/?"- (BEDT-TTF) 2 SF 5 CH 2 CF2S03 is a metal close to 
a metal-insulator transition. This assertion is corrob- 
orated by the fact that an external perturbation, such 
as a magnetic field |53| or pressure j34|, or replacing the 
anion CH2CF2 by CH2 pq |, drives the system into an 
insulating phase. 

Optical conductivity measurements on the 0-type crys- 
tals all show a broad feature from around 1000 to 5000 
cm" 1 H [TJ, [[J, EH- We assi g n tnis to the "mid- 
infrared" band. We now compare this interpretation to 
previous work. 6»-(BEDT-TTF) 2 CsZn(SCN) 4 , is a metal 
above 20 K and an insulator below 20 K. The frequency- 
dependent conductivity from 650 to 5000 cm -1 at tem- 
peratures above 20 K shows a feature from 650 to about 
1200 cm -1 which Tajima et al. assigned to a Drude 
peak, with a phcnomcnological damping rate of 0.1 eV 
~ 800 cm- 1 . 6>-(BEDT-TTF) 2 I 3 is a metal which under- 
goes a transition to a superconductor at 3.6 K. Tamura 
et al. ]3(| measured the optical conductivity from about 
700 to 4000 cm -1 in the metallic phase. They assign the 
spectrum below about 1000 cm- 1 to the tail of a Drude 
peak with a phenomenological damping rate of 400 to 800 
cm -1 , depending on the temperature. We disagree with 
the assignment of these low energy features to a Drude 
peak because it requires large scattering rates, compa- 
rable to t, implying a "bad" metal. In the k materials 
the observed widths of the Drude peak (and thus the 
scattering rate) at low temperatures are orders of mag- 
nitude smaller than this, typically of the order of tens 
of cm- 1 |37|l . We would assign all of the spectrum ob- 
served for 6»-(BEDT-TTF) 2 I 3 || to the "mid-infrared" 
band. The experiment of Tamura et a/.|3^1 did not go 
to low enough frequencies to observe the actual Drude 
peak. The feature in (9-(BEDT-TTF) 2 CsZn(SCN) 4 that 
Tajima et al. fl2]| assigned to a Drude peak we would 
assign to a phonon (and possibly its interference with 
the tail of the "mid-infrared" band) for the following 
reason. The optical conductivity spectrum of #-(BEDT- 
TTF) 2 RbCo(SCN) 4 was studied in detail by Tajima et 
al. |1J] . This material becomes a charge ordered insulator 
below 190 K. In slowly cooled samples the only signifi- 
cant spectral weight below 2000 cm -1 is a phonon around 
1200 cm- 1 (compare Fig. 5(a) in Ref. [TJ]). If the sample 
is quenched (i.e., rapidly cooled) there is a broad feature 
around this frequency. In the organics such quenching 
is usually associated with significant amounts of disorder 




Tajima et al. [jl2j assigned the second broad feature 
extending from above about 1200 cm -1 to domains asso- 
ciated with the onset of charge ordering. This was based 
on the calculated optical conductivity for a mean-field 
solution of the extended Hubbard model. Although we 
have essentially the same assignment as theirs for the 
" mid-infrared" band we note that their mean-field calcu- 



lation cannot actually produce this feature in the metallic 
phase. This underscores the many-body physics underly- 
ing this feature and the need for our exact diagonalization 
calculation. 

We suggest two possible reasons for the absence of the 
Hubbard resonance in the experimental spectra. First, 
Fig. 3 (c) and (d) suggest that most of the spectral weight 
from it is transfered to the " mid- infrared" band Second, 
it could be at a frequency above the range of the experi- 
ments. The most reliable estimates put U — V ~ 0.5 eV 
~ 4000 cm- 1 ffJJ. If V ^ It ~ 0.1 eV fTJ, Fig. 3 (c) 
would put the mid-infrared band extending from about 
1000 cm- 1 to 4000 cm -1 and the Hubbard band feature 
roughly around 6000 cm -1 . The latter is near the edge 
of many of the experimental plots. 

In summary, we have shown that increasing the 
nearest-neighbors Coulomb repulsion, V, for the quarter- 
filled extended Hubbard model on the square lattice, 
leads to an insulating phase for V > 2.2t. The calculated 
optical conductivity spectra close to the metal-insulator 
transition is consistent with experimental data on several 
9 and fi"- type molecular crystals. At V « V,^ 11 , the 
Drude peak is strongly supressed and the spectral weight 
is dominated by a broad "mid-infrared" band associated 
with short-range charge ordering. 

We thank O. Gunnarsson for making available his ex- 
act diagonalization code, E. Koch, M. Capone, F. Becca, 
M. Dressel, and M. Brunner for very helpful discussions, 
and the Max-Planck-Forschungspreis for financial sup- 
port. J.M. was supported by a Marie Curie Fellowship of 
the European Community programme "Improving Hu- 
man Potential" under contract No. HPMF-CT-2000- 
00870 and R.H.M. by the Australian Research council. 



* calandra@lmcp.jussieu.fr; Present address: Laboratoire 
de Mineralogie-cristallographie, Universite Pierre et 
Marie curie, 4 Place Jussieu, 75252 Paris, France 
[1] M. Vojta et al. Phys. Rev. B 62, 6721 (2000), and refer- 
ences therein. 

[2] M. B. Salamon et al. Rev. Mod. Phys. 73, 583 (2001). 

[3] S.K. Park et al. Phys. Rev. B 58, 3717 (1998). 

[4] Y. Ueda et al, J. Alloys Comp. 317-318, 109 (2001); T. 

Yamauchi et al, submitted to Nature. 
[5] See for example, D.S. Chow et al., Phys. Rev. Lett. 85, 

1698 (2000). 

[6] For a recent review, see: T. Ishiguro et al. Organic Su- 
perconductors , Springer- Verlag, Berlin, 2nd Edn, 1998. 

[7] H. Mori et al, S. Tanaka, and T. Mori, Phys. Rev. B 57, 
12023 (1998); T. Imakubo et al, J. Mater. Chem. 12, 
159 (2002). 

[8] For the importance of V in quasi-one-dimensional organ- 
ics see J. Hubbard, Phys. Rev. B 17, 494 (1978), and S. 
Mazumdar and A. N. Bloch, Phys. Rev. Lett. 57, 207 
(1983), J. Favand and F. Mila, Phys. Rev. B 54, 10425 
(1996). 

[9] H. Seo, J. Phys. Soc. Jap. 69, 805 (2000). 



5 



[10] 



[11] 
[12] 

[13] 

[14] 
[15] 
[16] 

[17] 
[18] 



[19] 

[20] 

[21] 

[22] 

[23] 
[24] 



H. Kino and H. Fukuyama, J. Phys. Soc. Jap. 65, 2158 
(1996); R. H. McKenzie, Comments Cond. Mat. Phys. 
18, 309 (1998). 

J. Dong et al, Phys. Rev. B 60 , 4342 (1999). 

H. Tajima et al, Phys. Rev. B 62, 9378 (2000). 

N. L. Wang et al, J. Phys.: Condens. Matter 13, 5463 

(2001). 

J. Ouyang et al, Phys. Rev. B 63, 054301 (2001). 
R. H. McKenzie et al, Phys. Rev. B 64, 085109 (2001). 
J. Merino and R. H. McKenzie, Phys. Rev. Lett. 87, 
237002 (2001). 

E. Dagotto, Rev. Mod. Phys. 66, 763 (1994). 
For a discussion of the relevant parameter values for the 
actual materials see Ref. [l5|. The most direct empiri- 
cal determination of U — V is for the dimers of BEDT- 
TTF+ in the insulator (BEDT-TTF) 2 Mo 6 Oi 9 based on 
the charge transfer excitation seen in the optical conduc- 
tivity. G. Visentini et al., Phys. Rev. B 58, 9460 (1998), 
obtain U - V ~ 0.5 eV. 

S. Mazumdar, R. T. Clay, and D. K. Campbell, Phys. 
Rev. B 62, 13400 (2000). 

R. T. Clay, S . Mazumdar, and D. K. Campbell, cond 



mat/0112278 



W. Stephan and P. Horsch, Phys. Rev. B 42 8736 (1990); 
E. Dagotto et al., Phys. Rev. B 45, 10741 (1992). 
D. J. Scalapino, S. R. White, and S. Zhang, Phys. Rev. 
B 47, 7995 (1993). 

W. Kohn, Phys. Rev. 133, A171 (1964). 

We find that, unlike the Drude weight, the charge gap: 



[25] 
[26] 



A c = E (N+l)+E (N-l)-2E (N), is strongly affected 
by finite size effects. 

E. Brezin and J. Zinn- Justin, Nucl. Phys. B 257, 867 
(1985). 

This estimate possibly overestimates , since the 

Lanczos method is limited to L < 20 and a priori we 
cannot exclude large finite-size effects. Taking L — 16, 20 
in the extrapolation gives a lower estimate of the transi- 



[27] 

[28] 

[29] 
[30] 
[31] 
[32] 
[33] 
[34] 
[35] 
[36] 
[37] 



[38] 



tion V r c 



1.2t; compared to V c 



0.7 1, obtained by 



a calculation based on slave-bosons J15|] . 

M. Vojta, A. Hiibsch, and R. M. Noack, Phys. Rev. B 

63, 045105 (2001). 

Y. Ohta et al., Phys. Rev. B 50, 13594 (1994); Physica 
B 194-196, 1559 (1994). 

K. Penc and F. Mila, Phys. Rev. B 49, 9670 (1994). 

J. A. Schlueter et al, J. Mater. Chem. 11, 2008 (2001). 

J. Favand and F. Mila, Phys. Rev. B 54, 10425 (1996). 

P.F. Maldague, Phys. Rev. B 16, 2437 (1977). 

J. Wosnitza et al., Phys. Rev. Lett. 86, 508 (2001). 

J. Hagel and J. Wosnitza, private communication. 

B. H. Ward et al, Chem. Mater. 12, 343 (2000). 

M. Tamura et al., J. Phys. Soc. Japan 57, 3239 (1988). 

Kornelson et al., Solid State Com- 
M. Dressel et al, Phys. Rev. B 50, 



See for example, K. 
mun. 72, 475 (1989) 
13603 (1994). 
T. F. Stalcup, J. S. 



Brooks, and R. C. Haddon, Phys. 



Rev. B 60, 9309 (1999). 



